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Using gauge/gravity duality, we compute the energy density and angular distribution of the power 
radiated by a quark undergoing circular motion in strongly coupled N = 4 supersymmetric Yang- 
Mills (SYM) theory. We compare the strong coupling results to those at weak coupling, and find 
the same angular distribution of radiated power, up to an overall prefactor. In both regimes, the 
angular distribution is in fact similar to that of synchrotron radiation produced by an electron in 
circular motion in classical electrodynamics: the quark emits radiation in a narrow beam along its 
velocity vector with a characteristic opening angle a ~ I/7. To an observer far away from the 
quark, the emitted radiation appears as a short periodic burst, just like the light from a lighthouse 
does to a ship at sea. Our strong coupling results are valid for any strongly coupled conformal field 
theory with a dual classical gravity description. 



I. INTRODUCTION AND OUTLOOK 

The emission and propagation of radiation is an in- 
teresting and natural topic of study in any field theory. 
For localized sources, one is typically interested in the 
angular distribution and spectrum of radiation at spatial 
infinity and in deciphering how information about the 
source is encoded in it. This topic has been studied for 
more than a century within the realm of classical electro- 
dynamics (for a textbook treatment, see [T]). It is much 
less studied in realms where quantum effects can be im- 
portant, as for example in a nonabelian gauge theory in 
the limit of strong coupling. This has largely been due 
to the absence of strong coupling methods for studying 
real-time dynamics. With the discovery of gauge/gravity 
duality, however, the dynamics of quantum fields in cer- 
tain strongly coupled nonabelian gauge theories can be 
mapped onto the dynamics of classical fields living in 
a higher dimensional spacetime [2!. Because of this, 
strongly coupled real-time phenomena in theories with 
gravitational duals have become accessible to reliable cal- 
culation. 

To date, much of the analysis of real-time phenomena 
in gauge theories with dual gravitational descriptions has 
focused on systems at nonzero temperature. This is be- 
cause these strongly coupled plasmas are more similar 
to the strongly coupled quark-gluon plasma of QCD be- 
ing produced and studied in heavy ion collision experi- 
ments at the Relativistic Heavy Ion Collider j3Hlj than 
the vacua of these gauge theories are to the vacuum of 
QCD. Asking how a test quark moving along some tra- 
jectory with nonzero acceleration radiates in vacuum in 
a strongly coupled nonabelian gauge theory is a simpler 
question than many of those that have already been stud- 
ied at nonzero temperature. We shall answer this ques- 
tion for the specific case of a test quark in circular mo- 
tion in the vacuum of M — 4 supersymmetric Yang-Mills 
(SYM) theory. 

Af = 4 SYM theory is a conformal nonabelian gauge 
theory, specified by two parameters: the number of colors 
N c and the 't Hooft coupling X = g 2 N c with g 2 the gauge 



coupling. We shall work in the large N c limit through- 
out. Since the theory is conformal, we can specify A. If 
we choose A — > 0, we can study the radiation using meth- 
ods similar to those used in classical electrodynamics. If 
we choose A — > 00, the dual gravitational description of 
this strongly coupled nonabelian quantum field theory 
becomes classical. Specifically, it becomes the classical 
gravity limit of Type IIB string theory in the AdSs x S5 
spacetime. The geometry of the five-sphere will play no 
role in our computation — indeed it can be replaced by 
any compact five-dimensional space X 5 . This means that 
our strong coupling results will be valid for all confor- 
mal quantum field theories with a dual classical gravity 
description — since conformality of the quantum the- 
ory maps onto the presence of an AdSs spacetime in the 
gravitational description. 1 

Because the vacuum of TV = 4 SYM theory is not sim- 
ilar to that of our world, studying the radiation of a test 
quark in circular motion in this theory does not have 
direct phenomenological motivation. However, we shall 
find that the angular distribution of the radiation is very 
similar to that of classic synchrotron radiation: when 
the quark is moving along its circle with an increasingly 
relativistic velocity, the radiation is produced in an in- 
creasingly tightly collimated beam — a beam which, if 
Fourier analyzed, contains radiation at increasingly short 
wavelengths and high frequencies. This is interesting first 
in its own terms and second with a view to extending 
the calculation to nonzero temperature. First, from the 
point of view of the quantum field theory we do not un- 



Our results are valid for all strongly coupled conformal quantum 
field theories with a dual classical gravity description as long as 
we express our results in terms of the 't Hooft coupling A in 
each of the conformal field theories. The relations between this 
parameter, N c , and the parameters that specify the gravitational 
physics in the AdSs space — namely the string coupling and 
the dimensionless ratio of the square of the AdS curvature and 
the string tension — will be different in different conformal field 
theories, since these relations do depend on the geometry of X5. 
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derstand why the qualitative properties of the radiation 
are so similar at A — > oo and A — > 0. Why should the 
pattern of radiation when what is radiated is nonabelian 
and strongly coupled be similar to that seen in classi- 
cal electrodynamics? In this sense, our results come as 
something of a surprise to us. Second, our results in this 
formal setting nevertheless open the way to a new means 
of modelling jet quenching in heavy ion collisions. If we 
were to add a nonzero temperature to our calculation, we 
could watch the tightly collimated beam of synchrotron 
radiation interact with the strongly coupled plasma that 
would then be present. The beam of radiation should be 
slowed down from the speed of light to the speed of sound 
and should ultimately thermalize, and it would be pos- 
sible to study how the length- and time-scales for these 
processes depend on the wavelength and frequency of the 
beam. 

We shall consider a test quark moving on a circle of 
radius Rq with constant angular velocity lu . In spherical 
coordinates {r, 9, tp} we can take the quark's trajectory 
to be given by 



7T 

Rq, 0=2' V = Uvt 



(1.1) 



By a test quark we mean a test charge in the fundamen- 
tal representation of the SU(N C ) gauge group. The total 
power radiated by this test quark moving in the vacuum 
of Af = 4 SYM theory in the large- N c and strong cou- 
pling limit has been computed previously by Mikhailov 
in Ref. [7J, who finds 



p 

P = a 

2tt 



(1.2) 



where a is the quark's proper acceleration, given by 



a = 7 vujq 



where the quark has speed 



and where 
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RqUq 



(1.3) 



(1.4) 



(1.5) 



is the standard Lorentz factor. Mikhailov observed that 
upon making the substitution 



(1.6) 



the result (1.2) becomes identical to the Larmor result 
in classical electrodynamics. Furthermore, Mikhailov 
showed that the power radiated by a test quark moving 
along an arbitrary trajectory is identical to what Lienard 
calculated in 1898 j8] — see [T] for a te xtbook treatment 
— upon making the substitution (1.6). We shall go be- 



yond Mikhailov's result by computing the distribution in 



space and time of the power radiated by the quark in 
circular motion. We shall not attempt Mikhailov's gen- 
eralization to an arbitrary trajectory. 

Before leaping to the conclusion that Mikhailov's re- 
sult suggests that the Af = 4 SYM synchrotron radiation 
pattern will be similar to that in classical electrodynam- 
ics, consider the results in the literature that push one's 
intuition in the opposite direction. One interesting study 
related to zero temperature radiation was performed in 
Ref. [9] , where decays of off-shell bosons in strongly cou- 
pled holographic conformal field theories were analyzed. 
(One can think of the off-shell boson as the virtual pho- 
ton produced in electron-positron annihilation, with that 
photon coupling to matter in a strongly coupled confor- 
mal field theory instead of to quarks in QCD.) It was 
shown in Ref. [9] that there is no correlation between the 
boson's spin and the angular distribution of power radi- 
ated through the sphere at infinity. In other words, if one 
prepares a state containing an off-shell boson with defi- 
nite spin, in the rest frame of the boson the event aver- 
aged angular distribution of power at infinity will always 
be isotropic. This stands in stark contrast to similar weak 
coupling calculations in QCD, where the boson's spin im- 
prints a distinct "antenna" pattern in the distribution of 
power at infinity j!0H12j . Similar behavior regarding the 
isotropization of radiation at strong coupling was also 
reported in Refs. [Ml [14]. A natural question to ask 
when considering the results of Ref. [9] is where does the 
isotropization of radiation come from? Is isotropization 
a characteristic of the particular initial states studied in 
Ref. [9] or is it a natural process which happens during 
the propagation of strongly coupled radiation? Indeed, 
it has been suggested in Ref. [T3] that the mechanism re- 
sponsible for the isotropization is that parton branching 
is not suppressed at strong coupling and, correspond- 
ingly, successive branchings can scramble any initially 
preferred direction in the radiation as it propagates out 
to spatial infinity. Regardless of mechanism, these cal- 
culations provide an example in which the pattern of ra- 
diation is very different in the strongly coupled Af = 4 
SYM theory to that in weakly coupled QCD or QED. If 
the intuition about the emission and propagation of radi- 
ation derived from these calculations applies to radiation 
from an accelerated test charge, then, we should expect 
that in our strongly coupled nonabelian theory this radi- 
ation should isotropize and "branch" , populating longer 
and longer wavelength modes as it propagates. We find 
nothing like this. The straightforward guess based upon 
Mikhailov's result — namely that synchrotron radiation 
in strongly coupled Af = 4 SYM theory should look like 
that in classical electrodynamics — turns out to be close 
to correct. 

We shall analyze the angular distribution of power ra- 
diated by a rotating quark in Af = 4 SYM theory at both 
weak and strong coupling. Our weak coupling results are 
valid only in this theory. At strong coupling, we obtain an 
analytic expression for the energy density which is valid 
at all distances from the quark and for all holographic 
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FIG. 1: A cartoon of the gravitational description of syn- 
chrotron radiation at strong coupling. The arena where the 
gravitational dynamics takes place is the 5d geometry of 
AdSs. A string resides in the geometry and an endpoint of 
the string is attached to the Ad boundary of AdSs, which is 
where the dual quantum field theory lives. The trajectory of 
the endpoint of the string corresponds to the trajectory of the 
dual quark. Demanding that the endpoint rotates results in 
the string rotating and coiling around on itself as it extends in 
the AdSs radial direction. The presence of the string in turn 
perturbs the geometry and the near-boundary perturbation 
in the geometry induces a Ad stress tensor on the boundary. 
The induced stress has the interpretation as the expectation 
value of the stress tensor in the dual quantum field theory. 



conformal field theories. At weak coupling our analy- 
sis reduces to that of synchrotron radiation in classical 
electrodynamics with an additional contribution coming 
from a scalar field [T5]- At strong coupling, gauge/gravity 
duality maps the problem onto that of the linear gravi- 
tational response due to the presence of a string. 

Fig. [T] shows a cartoon of the classical 5d dual grav- 
itational description. The five dimensional geometry of 
AdSs, which is the arena where the dual gravitational 
dynamics takes place, has a four-dimensional boundary 
whose geometry is that of ordinary Minkowski space. 
Ending at the boundary is a classical string whose end- 
point follows a trajectory that corresponds to the trajec- 
tory of the quark in the dual field theory. The condition 
that the endpoint rotates about a specific (Minkowski 
space) axis results in the entire string rotating about the 
same axis with the string coiling around on itself over 
and over. In turn, the rotating motion of the string 
creates gravitational radiation which propagates up to 
the boundary of the geometry. Just like electromagnetic 
fields induce surface currents on conductors, the gravita- 
tional disturbance near the boundary induces a Ad stress 
tensor on the boundary [THl [17]. As shown in Fig. [T] 
the induced Ad stress tensor inherits the coiled structure 



of the rotating string. Via the standard gauge/gravity 
dictionary [2j [18] , the induced Ad stress tensor is the ex- 
pectation value of the stress tensor in the dual quantum 
field theory. So, by doing a classical gravitational calcu- 
lation in five dimensions (which happens to be somewhat 
analogous to a classical electromagnetic calculation) we 
can compute the pattern of radiation in the boundary 
quantum field theory, including all quantum effects and 
working in a strong coupling regime in which quantum 
effects can be expected to dominate. 

At both weak and strong coupling we find that the ra- 
diation pattern at infinity qualitatively resembles that of 
synchrotron radiation in classical electrodynamics. And, 
at both weak and strong coupling, for the case of rel- 
ativistic circular motion we find the characteristic pat- 
tern of a lighthouse-like beam propagating out to infinity 
without broadening. As stated above, at weak coupling 
the difference between our results and those of classi- 
cal electrodynamics simply amount to the presence of an 
additional contribution coming from a scalar field. Re- 
markably, up to a normalization, we find that both the 
total power and the time-averaged angular distribution 
of power are the same at weak coupling and strong cou- 
pling. 

As the radiation is far from isotropic at infinity, this 
problem provides a concrete example demonstrating that 
radiation need not isotropizc in a nonabelian gauge the- 
ory at strong coupling. And, the width of the outgoing 
pulse of radiation does not broaden as the pulse prop- 
agates outward, meaning that as it propagates out to 
infinity the radiation remains at the (short) wavelengths 
and (high) frequencies at which it was emitted. 

We close this introduction with a look ahead at ex- 
tending this calculation to nonzero temperature T. The 
authors of Rcf. [TH] have shown that the power that it 
takes to move the quark in a circle is given by the vac- 
uum result (1.2 1 even for a test quark moving through 



the strongly coupled plasma present at T ^ as long as 



w 2 7 3 » 7r 2 T 2 



(1.7) 



In the opposite regime, where u> 2 -f 3 <C ir 2 T 2 , the power 
it takes to move the quark in a circle through the plasma 
is the same as it would be to move the quark in a straight 
line with speed v. This power has been computed using 
gauge/gravity duality in Refs. [2DHH], and corresponds 
to pushing the quark against a drag force. The ana- 
logue of our calculation — namely following where the 
dissipated or "radiated" power goes — was done, again 
using gauge/gravity duality, in Refs. [2"2Tf2"5] . There, it 
was demonstrated that the test quark plowing through 
the strongly coupled plasma in a straight line excites 
hydrodynamic modes of the plasma, trailing a diffusion 
wake and — if the motion is supersonic — creating a 
Mach cone. It is a generic feature of a nonzero temper- 
ature plasma, at weak coupling or at strong coupling, 
that any power dumped into it in some localized fashion 
must eventually take the form of hydrodynamic excita- 
tions 1251 . These considerations and our results motivate 
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the extension of the present calculation to nonzero tern- N c limit the classical effective Lagrangian reads 



perature, in particular in the regime (1.7). At these high 



velocities or low temperatures, we know that (1.2) is sat 



isfied and, given our experience at zero temperature in 
this paper, we therefore expect that at length scales that 
are small compared to 1/T the rotating quark radiates 
a narrow beam of synchrotron radiation, even at strong 
coupling. If (1.7) is satisfied, the wavelengths charac- 



terizing the pulse of radiation heading outward into the 
plasma are much shorter than 1/T. But, we also know 
on general grounds that this narrow pulse of beamed ra- 
diation must first reach local thermal equilibrium, likely 
converting into an outgoing hydrodynamic wave moving 
at the speed of sound, and must then ultimately dissipate 
and thermalize completely. Watching these processes oc- 
cur may give insight into jet quenching in the strongly 
coupled quark-gluon plasma of QCD. 

The outline of the rest of our paper is simple: we 
do the weak coupling calculation in Section II and use 
gauge/gravity duality to do the strong coupling calcula- 
tion in Section III. In Section IV we discuss the results 
and conclude. 



^classical = —F^F^ -J- A- ^(d X ) 2 - PX , (2-1) 

where F^ v is the field strength corresponding to the 17(1) 
gauge field A 11 , X is the scalar field, and the sources cor- 
responding to the test quark are given by 

p = e eS y/l - v 2 8 3 (r - r quark ) , (2.2) 

JM = ecS ^dT s3{r ~ rquark) ' (2 - 3) 

with r-q Uark the trajectory of the test quark, v = r^ uark 
the velocity of the test quark, and where it turns out that 
both effective couplings take on the same value 



' 2N r 2 



(2.4) 



In classical electrodynamics, p = and, in the source 
for the vector field J M , e e g is just e. We shall take the 
trajectory of the test quark to be ( 1.1 1 throughout. 



II. WEAK COUPLING CALCULATION 



A. Solutions to the equations of motion and the 
angular distribution of power 



In the limit of weak coupling, the radiation produced 
by an accelerated test charge can be analyzed using clas- 
sical field theory. For example, for the case of an electron 
undergoing synchrotron motion in QED, the appropriate 
classical effective description is simply that of classical 
electrodynamics. For the case of a test quark (specifi- 
cally, an infinitely massive particle from the Af = 2 hy- 
permultiplet that is in the fundamental representation of 
the SU(N C ) gauge group) moving through the vacuum of 
JV = 4 SYM theory, more fields can be excited and the 
appropriate effective classical description is more com- 
plicated. The Lagrangian for this theory can be found, 
for example, in Ref. [26] . As can readily be verified by 
inspecting the field theory Lagrangian, in the limit of 
asymptotically weak coupling and large test quark mass 
the fields that can be excited by the test quark are the 
non-Abelian gauge field and an adjoint representation 
scalar field. Interactions with all other fields are either 
suppressed by an inverse power of the quark mass or a 
positive power of the 't Hooft coupling A. Moreover, in 
the limit of asymptotically weak coupling the gauge and 
scalar fields satisfy decoupled linear equations of motion. 

Following Ref. fT5], instead of dealing with N% decou- 
pled components of the gauge field and N% decoupled 
components of the adjoint scalar field, we consider the 
case of a single U(l) gauge field and a single real scalar 
field with appropriate effective couplings tailored to com- 
pensate for the fact that the actual theory contains ad- 
joint representation fields. The effective couplings can be 
determined by computing the vector and scalar contribu- 
tions to the Coulomb field of a quark in the effective de- 
scription and matching onto the full theory. In the large 



Comprehensive analyses of the radiation coming from 
the gauge field can be found in classical electrodynamics 
textbooks (see for example p]). However, because the 
case of scalar radiation is less well known and also for the 
sake of completeness, we present a brief analysis of the 
radiation for both fields. We shall linger in our discussion 
on those aspects which will also be of value in the analysis 
of our strong coupling results. 

Choosing the Lorentz gauge d ■ A = 0, the equations 
of motion for the fields read 



- d 2 M l = , -d\ = P ■ 
These equations are solved by 



(2.5) 



A"(t,r) = d i r , G(t-t',r-r , )J f *{t , ,r r ), (2.6a) 



X(t,r) = 

where 

G{t~t',r-r') 



d A r'G(t~t',r-r')p{t',r'), (2.6b) 



6{t-t') 



4tt\i — r 



' 6(t-t'-\r-r>\) (2.7) 



is the retarded Greens function of — d 2 , and (f , r') is the 
point of emission. After using delta functions to carry 
out the r' M integrations, we arrive at the expressions 



X 



e e ff 1 Rquark 

47rrS(i',r) dt' 
e e ff 1 1 



t'=tret 



4nr E(i', r) 7 



t'=*re 



(2.8a) 
(2.8b) 
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FIG. 2: A sketch of the solutions (2.81. As the quark moves 



along its circular trajectory, it emits radiation in a narrow 
cone of angular width a ~ I/7 in the direction of its velocity 
vector. The diagram is a snapshot at the time when the quark 
is at the top of the circle. The red spiral shows where the 
radiation emitted at earlier times is located at the time of the 
snapshot. The width A of the spiral scales like A ~ I/7 3 , as 
explained in Fig. [3] 



where 



"(£' r) = l r r quark(^ )| r ' ^quarkt/ ) ^ 



7 = — 1> 2 , r = \r\, and the retarded time i re t is the 

solution to 



t - trot — \r — r quark (i rot )| = . 



(2.10) 



(t—t Ic t) is the light travel time between r quar k(iret) and r. 
Because the motion of the quark is periodic, r qU ark(*ret) 
is a periodic function of its argument with period 2tt/ujq. 
Eq. (2.10) then implies that if t is shifted by 2tt/ujq, t lct 
shifts by the same amount so that (t — i re t) is left un- 
changed. Because the motion of the quark is circular, 
Eq. (2.101 also implies that a change in the position of 



the quark corresponding to a shift in its azimuthal angle 
by dip is equivalent to a shift in both t and t let by 8ip/u)Q. 

We now describe the qualitative behavior of the solu- 
tions (2.8). Fig. [2] shows a pictorial representation of the 
solutions ( |2.8[ ) . As the quark moves along its circular tra- 
jectory, it emits radiation in the direction of its velocity 
in a narrow cone of angular width a. From the solutions 
(2.8) it can be shown that in the large 7 limit a scales 
like a ~ I/7 for both the scalar and vector radiation. 
The cone of radiation emitted at each time propagates 
outwards at the speed of light. At any one moment in 
time, the radiation emitted at all times in the past forms 
a spiral, as illustrated by the red spiral sketched in Fig. [2] 
Clearly, the radial width A of the spiral must go to zero 
as a — > 0, namely as 7 — > 00. As we illustrate in Fig. [3j 
in the large 7 limit the radial width scales like A ~ 1/ 7 . 
To understand the scaling of A, consider an observer at 
the point p in Fig. [3j As the quark moves along its tra- 
jectory, an observer at p will see a short pulse of radiation 




FIG. 3: A close-up illustration of the emission of radiation at 
two times ti and t^. The radiation is emitted in the direction 
of the quark's velocity vector, within a cone of angular width 
a. An observer at the point p is illuminated by a pulse of 
radiation of duration At ~ Roa and of spatial thickness A. 
The leading edge of the pulse observed at p is emitted at ti 
and the trailing edge observed at p is emitted at t% . At time £2 
the radiation emitted at t\, denoted by the solid red line, has 
traveled a distance Rou/v towards p. The chordal distance 
between the two emission points is Roa in the a — > limit. 
The width A is therefore A = R a(l/v - 1). 



of duration At ~ R^a/v. The leading edge of the pulse 
will be emitted by the quark at time t\ and the trail- 
ing edge will be emitted at time ti = t± + R^a/v. At 
time ti the radiation emitted at t\ will have traveled a 
distance R a/v towards p. Moreover, the chordal dis- 
tance between the two emission points is approximately 
Roa when a is small. It therefore follows that the spatial 
thickness A of the pulse observed at p scales like 



(R a/v - R a a) - i? a/ 7 2 - i? /7 3 



(2.11) 



We note that the fact that A and a scale with different 
powers of 7 is a consequence of the fact that the radi- 
ation is emitted in the direction of the quark's velocity. 
If the radiation were emitted in any other direction, A 
and a would have the same 7 scaling. Because we will 
re-use these arguments in the analysis of our strong cou- 
pling results, it is important to note that they are purely 
geometrical, relying only on the fact that the radiation 
travels at the speed of light and the fact that the outward 
going pulse of radiation propagates without broadening, 
maintaining the spiral shape whose origin we have illus- 
trated. 

In addition to determining how the pulse widths a 
and A are related, these geometrical considerations in 
fact specify the location of the spiral of radiation pre- 
cisely. For a quark moving in a circle of radius Rq in the 
(x, y) plane, at the time at which the quark is located at 
(x, y) = (Rq,0) the spiral of energy density radiated at 
all times in the past is centered on a curve in the z = 



G 



plane specified by 



x = Rq 



cos p+ — sm p 

v 



sm p+— cosp 
v 



y = Ro 

r 

(2.12) 

where p is a parameter running from to oo. At large r, 
this spiral can equivalently be written 



(p (r,t) = -+uj t- — +O(l/r) 
I Ho 



(2.13) 



in spherical coordinates, where we have also added the 
i-dependence. 

Recall that the 7-scaling of A depended on the fact 
that the radiation is emitted tangentially, in the direc- 
tion of the velocity vector of the quark. The same is true 
for the shape of the spiral. For example, if we repeat the 
geometrical analysis for a hypothetical case in which the 
quark emits a narrow cone of radiation in the f direc- 
tion, perpendicular to its direction of motion, the energy 
density would be centered on a different spiral specified 
by 



X = Rq 



1 



9 



cosp, y = -Rq 



1 



sin p. (2.14) 



This spiral, including its i-dependence, is given in spher- 
ical coordinates by 



(Technically, the integration in Eq. (2.16) should run 
from t Ie t — c to t Ie t = c + 2tt/ujq where c ~ — r. This 
is simply because in the far-zone limit the retarded time 
always scales like t YCt ~ — r. However, as the motion of 
the quark is periodic one can always shift the integration 
variable so that the integration runs from to 2tt/luq.) 

For the simple case of scalar or electromagnetic fields, 
the computation of the energy flux directly from the fields 
is relatively simple. However, in the gravitational set- 
ting discussed below, the computation of the energy flux 
requires solving three differential equations whereas the 
computation of the energy density requires only solving 
one. Therefore, if possible, it is useful to extract the an- 
gular distribution of power from the energy density alone. 

For a conformal theory at zero temperature the far- 
zone energy density and flux associated with a localized 
source must take the form 

£(t,r) = ^M S(t,r)= S ^ 2 6 ^ r, (2.18) 

up to 0(l/r 3 ) corrections. Moreover, using the continu- 
ity equation and the fact that dt le t/dt = —dt re t/dr in 
the far zone, it is easy to see that we must have e = s. 
In terms of e, the time-averaged angular distribution of 
power reads 



ip(r, t) = v + uj t 



Ra 



(2.15) 



dP oj r <"o dt 
dQ 2tt J q rct dtret 



(2.19) 



which differs from the correct spiral of Eqs. (|2.12|) and We see from Eq. ( |2.19[ ) that we can extract the angular 



( 2.13 ) by a phase shift at large r. If the quark emits radi- 



ation perpendicular to its direction of motion, a distant 
point p would be illuminated by light emitted when the 
quark was at the closest point on its circle to p, whereas 
for the correct spiral (2.12) the light seen at p was emit- 
ted one quarter of a period earlier but had travelled Rq 
farther, leading to the phase shift of (— tt/2 + v). In the 



distribution of the radiated power at infinity by analyzing 
the asymptotic behavior of the energy density. 

The energy density associated with the vector field is 
^vector = \{E 2 + B 2 ) where E and B are the electric and 
magnetic fields. The energy density associated with the 
scalar is £ sca iar = \ \{dtX) 2 + ( ^x) 2 ] • Evaluating these 
quantities for the solutions ( 2.8 1 in the r — > 00 limit yields 



near zone, ( 2.14[ ) differs from the correct spiral (2.12) in 
shape, not just by a phase shift. 

In Section IV, we shall compare our results at strong 
coupling to those we have obtained in this Section via ge- 
ometrical arguments built upon the fact that the outward 
going pulse of radiation propagates without spreading. 

We now analyze the far zone behavior of the fields 
and compute the time-averaged angular distribution of 
power. In particular, we wish to compute the angular dis- 
tribution of power radiated through the sphere at r = 00 
averaged over one cycle of the quark's motion. Denoting 
the energy flux by S, the time averaged angular distri- 
bution of power is given by 



3 + 2« 2 + (l-2v 2 )cos26» 



647r 2 r 2£6 

— 8v sin 9 sin(iy9 — wo^ret ) 

— 2 sin 2 9 cos (2(ip — wot re t)) 



(2.20a) 



-'scala 



^2 „2. ,2 



167r 2 7 2 r 2£6 



sin 2 6»cos 2 ((^-w trct), (2.20b) 



up to subleading corrections proportional to 1 /r 3 , where 
we have defined 



dP 



dfl 2?r r- 



U y 2 

lim r 



dt r 



dt 

dt T ei 



where, from Eqs. (2.10) and (2.9), 



dt 



9t rR t t — t T 



(2.16) 



(2.17) 



£ = 1 — v sin 9 sin(</? — woiret) 
which is the r —¥ 00 limit of S: 

5 = £ + 0(l/r). 



(2.21) 



(2.22) 



Note that the energy densities (2.20) are proportional to 
l/£ 6 and therefore have maxima in the 9 — tt/2 plane 



7 



where £ has minima; at large r, these occur on the spiral 

Noting from (2.10) that (f — i rct ) ~ r in the r — » oo 
limit, we see from (2.17) and (2.22) that dt/dt lct = £ 



in the r — > oo limit. We can now evaluate the expres- 
sion (2.19) for the time-averaged angular distribution of 



power, yielding 



dP vl 



2 9 2 



8-4 sin 2 6» - v 2 (l + 3v 2 ) sin 4 6> 



scala 



1287T 2 



w 2 w 2 e 2 ff sin 2 i 



<; 2 sin 2 0) V2 



(4 + u 2 sin 2 9) 



dfl 



128tt 2 7 2 



2 c:„2 /A 7/2 



w z sm' 



(2.23a) 
(2.23b) 



Eqs. ( 2.23 1 show that the power radiated in vectors and in 



scalars are both focused about 9 = tt/2 with a character- 
istic angular width 59 oc 1/7. (Since 1 — v 2 sin 2 9 = I/7 2 
at 6 = tt/2 and l~v 2 sin 2 9 ~ (l + c)/7 2 at 9 = vr/2-c/7 
for any constant c when 7 > 1.) Adding the contribu- 
tions (2.23a) and (2.23b), we find that the time-averaged 



angular distribution of power radiated in vectors plus 
scalars is 



dPtotai _ « 2 Woe 2 ff 



9 • 2 

v sin 



9) 



327T 2 



[1 - v 2 sin 



.5/2 



(2.24) 



Integrating Eqs. (2.23) over all solid angles we find the 



power radiated in vectors and scalars 



P vector 
-^scalar 



2 p 2 

_feff 2 

3 4tt 
1 p 2 

3 4tt 



(2.25a) 
(2.25b) 



where a = vujq/ (1 — v 2 ) is the proper acceleration of the 



test quark. With e c s = e, Eq. (2.25a) becomes the well 



known Larmor formula for the power radiated by an ac- 
celerated charge in classical electrodynamics. Evidently, 
at weak coupling two times as much energy is radiated 
into the vector mode as the scalar mode. We now see 



that Mikhailov's strong coupling result (1.2) is equiva- 
lent to the result for weakly coupled N = 4 SYM theory 
upon making the substitution 



2 

ciT 



(2.26) 



with the difference between this result and ( 1.6 ) — which 
was obtained by comparing the strong coupling N = 4 
SYM result to classical electrodynamics — coming from 
the fact that in weakly coupled J\f = 4 SYM theory the 
rotating test quark radiates scalars as well as vectors. 



III. GRAVITATIONAL DESCRIPTION AND 
STRONG COUPLING CALCULATION 

Gauge/gravity duality [5] provides a description of 
the (strongly coupled; quantum mechanical) vacuum of 



J\f = 4 SYM theory in N c —> 00 and A —> 00 limit in 
terms of classical gravity in the AdSs x S5 spacetime. 
As we discussed in Section I, generalizing from J\f = 4 
SYM to any other conformal quantum field theory that 
has a classical gravity dual involves replacing the S5 by 
some other compact five-dimensional space. In what fol- 
lows, we shall never need to specify this five-dimensional 
space, meaning that all our strong coupling results will 
be valid for any large-7V c conformal field theories with 
gravitational duals. 

The metric of the AdSs spacetime is 



ds 2 = 



L 2 



-fdt 2 + dr 2 



dv 2 ^ 



(3.1) 



where L is the AdS$ curvature radius, u is the (inverse) 
AdS§ radial coordinate with the boundary at u — 0, and 
/ = 1. (We have kept the factor of / in the metric in 
order to generalize some of our calculations to nonzero 
temperature, where 



/ = 1 



(3.2) 



with Uh = 1/(tvT) and T the temperature of the field 
theory state.) A single test quark in the fundamental 
representation in the field theory corresponds in the grav- 
itational description to an open string with one endpoint 
attached to the boundary at u = 0. Taking the iV c — > 00 
and A — > 00 limits eliminates fluctuations of the string 
and makes the probability for loops to break off the string 
vanish. Upon taking these limits, the string is therefore 
a classical object. The presence of the classical string 
hanging "down" into the five-dimensional bulk geometry 
perturbs that geometry via Einstein's equations, and the 
behavior of the metric perturbation near the boundary 
encodes the change in the boundary stress tensor. This 
is illustrated in Fig. [T] and in this Section we shall turn 
this cartoon into calculation. 



A. The rotating string 

We begin by describing the rotating string. The metric 
perturbations that this string creates are small in the 
large N c limit, small enough that we can neglect their 
back reaction on the shape of the string itself. We can 
therefore begin by finding the shape of the rotating string 
in AdSs, in the absence of any metric perturbations. The 
dynamics of classical strings are governed by the Nambu- 
Goto action 



'NG 



= -Tn 



(3.3) 



where To = \/A/27rL 2 is the string tension, a and r 
are the coordinates on the worldsheet of the string, and 
g = det g a i, where g a j, is the induced worldsheet metric. 
The string profile is determined by a set of embedding 



functions X M (r, a) that specify where in the spacetime 
described by the metric Gmn the point (r, a) on the 
string worldsheet is located. The induced world sheet 
metric is given in terms of these functions by 



g ab = d a X ■ d b X . 



(3.4) 



where a and b each run over (t, a). For the determinant 
we obtain 



- g = (d T x ■ d a xf - (d T x) 2 (d„xy 



(3.5) 



We choose worldsheet coordinates r — t and a = u. 
As we are interested in quarks which rotate at constant 
frequency luq about the z axis, we parametrize the string 
embedding functions via 



X M (t,u) = (t,r s (t,u),u) 



(3.6) 



where in spherical coordinates {r, 9, if} the three- vector 
r s is given by 



,(t,«) = (R(i 



w 



(3.7) 



The function <j)(u) describes how the purple string in 
Fig. [T] winds in azimuthal angle as a function of "depth" 
in the AdSs radial direction, whose coordinate is u. The 
function R(u) describes how the string in Fig. [I] spreads 
outward as a function of depth. Because we are describ- 
ing a situation in which the test quark at the boundary 
has been rotating on its circle for all time, the purple 
string in Fig. [T] is rotating at all depths with the same 
angular frequency w , meaning that the two functions 
R(u) and 4>(u) in the parametrization (3.7) fully specify 
the shape of the rotating string at all times. With this 
parameterization, the Nambu-Goto action reads 



S 



NG 



VA 

2^r 



dtduC, 



where 



v /(l-^i? 2 )(l + .R' 2 )+i?V : 



(3.8) 



(3.9) 



To determine the shape of the string in Fig. [TJ we must 
extremize (3.8). 

(u) and R(u) follow from 
One constant 



The equations of motion for 
extremizing the Nambu-Goto action (|3.8[) 



of motion can be obtained by noting that the action is 
independent of (f>(u). This implies that 



n : 



dC 



(3.10) 



is constant. Here and throughout, by ' we mean d/du. 
The minus sign on the right-hand side of (3.10) is there 
in order to make II positive for positive loq. 
II, the equation of motion for '(u) reads 



In terms of 



u 4 n 2 (i-w 2 i? 2 )(i + i?' 



R 2 (R 2 



u 4 n 2 



(3.11) 



= 0. 



(3.12) 



The equation of motion for R(u) is given by 

dC _ d dC 
OR " d^~dR' 

Taking the above partial derivatives of C and then elimi- 
nating </> derivatives via Eq. (3.11), we obtain the follow- 
ing equation of motion: 

R{u + 2RR')(l + R' 2 ) l + R' 2 



R" 



i(u 4 n 2 - r 2 ) 



i?(l 



u 2 R 2 ) 



0. 

(3.13) 



As a second order differential equation, one might ex- 
pect that Eq. (3.13) requires two initial conditions for 



the specification of a unique solution. However, it was 
shown in Ref. UHl that this is not the case — together, 



Eqs. (3.11 ) and (3.13) imply that all required initial data 
for (3.13) is fixed in terms of II and ljq. To see how 



this comes about note that Eq. (3.13) is singular when 
or when R 2 — u 4 H 2 — 0. In order to main- 



1 



LO Q H 



tain the positivity of the right hand side of Eq. (3.11), 



uj 2 R 2 and R 2 



1 



This happens where 



1 II 2 must vanish at the same point 



Vfiwo' 



R — R c 



1 

u ' 



(3.14) 



We can then solve (3.13) in the vicinity of u = u c by 
expanding R(u) as 



1 



R{u) = R C + R' c {u - u c ) + -JZ£(u - u c ) 2 + 



(3.15) 



where R, 



O) 



d"R\ 



. The second and third terms 
in (3.13) are divergent at u — u c while the R" term is 



finite there. This means that when we substitute (3.15) 
into (3.13) and collect powers of (u — u c ), we obtain an 
equation for R' c that does not involve i?", meaning that 
the equation of motion itself determines R' c . A short 
calculation yields 



R' 



UJQ 

n 



(3.16) 



Therefore, both R c and R' c are determined by luq and 
II and consequently, for given ujq and II there is only 
one solution to Eq. (3.13). From the perspective of the 



boundary field theory this had to happen: for a given 
angular frequency wo, the radiation produced by the ro- 
tating quark can only depend on the radius Rq of the 
quark's trajectory (which as we will see below can be 
related to II). 

W ith th e above point in mind, the unique solution to 
Eq. (|3.13|) is given by 2 



R(u) 



V 2 "/ 2 U 2 



R i 



(3.17) 



One way to derive the solution | |3.17[ l is simply to solve Eq. l |3.13| l 
with the series expansion ( |3.16| l and recognize the resulting series 
as that of 13.171. 
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where v = Rqujq is the velocity of the quark and 7 
1 / \fl — v 2 . The velocity of the quark is related to II via 



TT 2 4 

11 = V 7 UJq ■ 



(3.18) 



With R(u) given by (3.17), it is easy to integrate 



Eq. (3.11). Taking the negative root of Eq. (3.11) so 



that the string trails its u = endpoint, we find 

4>{u) = —ujluo + arctan (ujujo) (3.19) 

and thus 

ip(t, u) = tuo(t — 7U) + arctan (wjujq) . (3.20) 

We thus have an analytic description of the shape of the 
rotating string. 

We can obtain the profile of the rotating str ing p ro- 
jected onto the plane at the boundary by using (3.17) to 
eliminate u from (3.20), yielding 



<p(t,r) = uj t- 



7? 2 



1 + arctan 




In the limit of r Rq , we have 
I K 



(3.21) 



(3.22) 



Equation (3.22 1 describes a spiral whose neighboring cir- 
cles (at fixed t) are separated by 2ttRq. Note that the 



curve ( |3.22 ) is the same spiral as (2.13) when the veloc- 
ity of the quark v 1. Indeed, it can also be shown that 
at all radii (3.21 ) is the same spiral as (2.12) with v = 1. 



We conclude that for v — > 1 the rotating string hanging 
down into the AdSs geometry is lined up precisely below 
the spiral in the boundary energy density in the limit of 
weak coupling and, as we shall see below, in the limit of 
strong coupling. 3 At velocities less than 1, the shape of 
the rotating string is not related to the shape of the spiral 
of energy density, at least not directly by projection. This 
is perhaps not surprising, as it is only in the v — > 1 limit 
that the widths (in radius, azimuthal angle, and polar 
angle) of the spiral of radiation emitted by the rotating 
quark all go to zero and the radiation becomes localized 
on the spiral (2.12). For v < 1, the spiral of radiation 



3 This identification suggests interesting possibilities at nonzero 
temperature. Whereas at T = the rotating string has a profile 
R(u) such that it extends out to infinite r, when T ^ the rotat- 
ing string only extends out to some finite maximum radius r m ax 
that depends on u>o and v 1191 . The fact that at zero temperature 
we can identify the the spiral of the rotating string with the spiral 
of energy density in the boundary quantum field theory suggests 
that, at least at V — y 1, the maximum radius out to which the 
rotating string reaches may be some length scale related to how 
far the beam of synchrotron radiation can travel through the hot 
strongly coupled plasma before it thermalizes. Interestingly, if 
one takes the v — > 1 at fixed Rq, one finds r max — > 00 I19| . 



is broadened and its location as a one-dimensional curve 
is no longer uniquely defined. In the v — >■ 1 limit, how- 
ever, the spiral of radiation and the spiral of the string 
in the bulk are both one-dimensional curves, and there 
must be some relation between these curves — a relation 
that turns out to be simple projection. 

The shape of the rotating string determines the total 
power radiated, as we now explain. The string has an 
energy density ir® and energy flux 7Tq, where 



5S- 



NG 



' M 



with <Sng given in (3.3), which satisfy 



d t TT° + CU* = 



«"0 



(3.23) 



(3.24) 



This is simply the statement of energy conservation on 
the string worldsheet. The energy flux down the rotating 
string is given by 



2tt 



w n 



2^r~w! 



2T 



(3.25) 



where we have used (3.18) and where a = vloqJ 2 is the 



proper acceleration of the test quark, and where the mi- 
nus sign comes from our choice of parametrization of the 
world sheet with a coordinate u that increases down- 
ward, away from the boundary. The energy flux down 
the string is the energy that must be supplied by the ex- 
ternal agent which moves the test quark in a circle; by 
energy conservation, therefore, it is the same as the total 
power radiated by the quark. 4 Thus far, we have repro- 
duced Mikhailov's result ( 1.2 ) via a calculation along the 
lines of Ref. [TH] — but with the advantage that we have 
been able to find analytic expressions for R(u) and </>(m), 
which will be a significant help below. 

We close this subsection by commenting on an aspect 
of the physics of the rotating string that we shall not 
pursue in detail. As pointed out earlier, there is a special 



The energy loss rate l |3.25| l can also be derived from the energy 
density along the string 

^ (1 + nuW 2 ) , (3.26) 
this energy density approaches the con- 



-0° 



as follows. For «»« c 
stant 



2n 



'ylluJO 



(3.27) 



Now, consider the total energy along a segment of the string 
whose extent in u is Su, which we shall denote by 6E = — tt^Su. 
The rotating quark can be thought of as "producing this much 
new string" in a time St that, from (3.20 1, is given by St = ^ . 
One should then identify SE with the energy pumped into the 



system during St. Equating SE ■■ 



-St = PSt we thus find that 



"To 



2tt 



IIa)o 



(3.28) 
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value of u, u = u c where we should specify the initial 
condition needed to solve the equations of motion that 
determine the shape of the rotating string. u c has a more 
physical interpretation: it is the depth at which the local 
velocity of the rotating string, ujqR(u c ), becomes equal 
to the speed of light in the five-dimensional spacetime. 
This has a simple interpretation from the point of view of 
the string worldsheet. It can be verifed that the induced 
metric on the string worldsheet g a b has an event hori- 
zon at u — u c . Thus disturbances on the string become 
causally disconnected across u c . (This suggests that in 
the boundary field theory the regime r 3> R(u c ) — Rq/v 
can be thought of as the far field region, while r < Ro/v 
can be thought of as the near field region.) Even though 
the AdSs spacetime with metric Gm n has zero tempera- 
ture, the worldsheet metric is that of a (l+l)-dimensional 
Schwarzschild black hole whose horizon is at u = u c and 
whose Hawking temperature is given by 



1 



2ttju c 



1, 

1 



Unruh 



(3.29) 



where Tunruh = is the Unruh temperature for an ac- 
celerated particle with proper acceleration a. As it ra- 
diates, the rotating test quark should experience small 
kicks which would lead to Brownian motion in coordi- 
nate space if the quark had finite mass. At strong cou- 
pling such fluctuations can be found from small fluctua- 
tions of the worldsheet fields X M (t, u) in the worldsheet 
black hole geometry. The thermal nature of the world- 
sheet metric indicates that such Brownian motion can be 
described as if due to the presence of a thermal medium 
with a temperature given by (3.29|. 



B. Gravitational perturbation set-up 

In the limit N c — > oo, the 5d gravitational constant 
is parametrically small and consequently the presence of 
the string acts as a small perturbation on the geometry. 
To obtain leading order results in 1/N C we write the full 
metric as 



G 



MN 



(0) 
MN 



IMN ■ 



(3.30) 



where G^l, is the unperturbed metric (|3.l|), and lin- 



earize the resulting Einstein equations in the perturba- 
tion Hmn ■ This results in the linearized equation of mo- 
tion 



-D 2 h 



MN 



2D P D( M h N ^p — D M D N / 



L 2 



L MN 



(D 2 h - D P D Q h PQ -£h)G 



— 2k\ Imn , 



MN 



(3.31) 



where h 



h M 

11 Mi 



Dm is the covariant derivative under 
the background metric (3.1 1, k 2 is the 5d gravitational 
constant and £mat is the 5d stress tensor of the string. In 
Af = 4 SYM theory, k\ = 4n 2 L 3 /N 2 , but this relation 
would be different in other strongly coupled conformal 
field theories with dual gravitational descriptions. We 
shall see that k 2 does not appear in any of our results. 

According to the gauge/gravity dictionary, the on-shell 
gravitational action 



Sc 




J2) + Sgh 



(3.32) 



is the generating functional for the boundary stress ten- 
sor [TTJ [TH] - Here, G is the determinant of Gmn, 71 is its 
Ricci scalar, and Sgh is the Gibbons-Hawking boundary 
term [57] , discussed and evaluated in the present context 
in Ref. [25 . The 5d metric Gmn induces a Ad metric 
on the boundary of the the 5d geometry. The boundary 
metric is related to the bulk metric by 



lim -prGav{x,u) 

u— >0 L 



(3.33) 



Because Gmn oc l/u 2 , the rescaling by u 2 in (3.33) yields 
a boundary metric that is regular at u — 0. The bound- 
ary stress tensor is then given by [17] 



2 SS G 



(3.34) 



with g denoting the determinant of g^. We see that in 
order to compute the boundary stress tensor, one needs 
to find the string profile dual to the rotating quark and 
compute its 5d stress tensor tMN- One then solves the 



linearized Einstein equations (3.31) in the presence of 



the string source and then extracts the boundary stress 
tensor fro m the v ariation of the on-shell gravitational ac- 



tion, as in (3.34). The Kg, dependence drops out because 



Sg oc l//t| while we see from (3.31 1 that the perturbation 
to the metric is proportional toKg. 

We determined the string profile in Section III. A. From 
this, we may now compute the 5d string stress tensor, 
which is what we need in order to determine the metric 
perturbation due to the string. In general, 



.MN 



gg^daX^'dbX'^S^r-rs). (3.35) 



For the rotating string given by (3.6) and (3.7| with 
(13. 171) and (13.191) , this reduces to 
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tMN — 



-J 3 (r - r s 



( R' 2 + R 2 s 2 (f>' 2 + 1 R 2 s 2 u Q R'(j)' -R 2 s 2 lu (R' 2 + 1) R 2 s 2 u cf>' \ 

R 2 s 2 Lo Q R'(f>' R ,2 (R 2 s 2 lu 2 -1) -R 2 s 2 R'fi R' (R 2 s 2 uo 2 - l) 



~R 2 s 2 u {R' 2 + 1) -R 2 s 2 R'<P' i? 4 s 4 (w 2 i?' 2 +w 2 -ft 2 ) -R 2 s 2 fi 

\ i? 2 sW R'(R 2 s 2 uj 2 - 1) -i? 2 sV i?Vw 2 -l / 



(3.36) 



r 



where the components are in spherical coordinates, with 
the rows and columns ordered (t,r,9,cp,u), and where 
s = sin#. 



Gauge invariants and the boundary energy 
density 



one can solve 
However, do- 



Given the string stress tensor ( 3.36[ ) 
the linearized Einstein equations ( 3.31[ ) 
ing this in its full glory is more work than is necessary. 
The metric perturbation Hmn contains fifteen degrees of 
freedom which couple to each other via the linearized Ein- 
stein equations. This should be contrasted with the Ad 
boundary stress tensor, which is traceless and conserved 
and thus contains five independent degrees of freedom. 
Therefore, not all of the degrees of freedom contained in 
h-MN are physical. The linearized Einstein equations are 
invariant under infinitesimal coordinate transformations 
X M X M + £ M where £ M is an arbitrary infinitesi- 
mal vector field. Under such transformations, the metric 
perturbation transforms as 



>>MN 



IAIN 



— Dm£n — -Dat£ 



M : 



(3.37) 



where Dm is the covariant derivative under the back- 
ground metric G A °j N . The physical information contained 
in Iimn must be gauge invariant. This requirement re- 
stricts the number of physical degrees of freedom in Hmn 
to five, 5 matching that of the boundary stress tensor. 

Just as in electromagnetism where gauge invariant 
quantities {e.g. electric and magnetic fields) can be con- 
structed from a gauge field, it is possible to construct 
gauge invariant quantities out of linear combinations of 
h,MN and its derivatives. The utility of doing this lies 



in the fact that the equations of motion for gauge in- 
variants don't carry any of the superfluous gauge vari- 
ant information contained in the full linearized Einstein 
equations and thus can be simpler to work with. Indeed, 
as demonstrated in Refs. [531 123 I2H]) the equations of 
motion for gauge invariants can be completely decoupled 
from each other. 

As the background geometry is translationally invari- 
ant, it is useful to introduce a Ad spacetime Fourier trans- 
form and work with mode amplitudes /imjv(w; w, q). Use- 
ful gauge invariants can then be constructed out of linear 
combinations of hMN{u;ui,q) and its radial derivatives 
and classified according to their behavior under rotations 
about the q-axis. As 1~imn is a spin two field, there ex- 
ists one independent helicity zero combination and a pair 
each of independent helicity one and two combinations. 6 

For utility in future work, we shall give an expression 
for the gauge invariant quantity that is valid at zero or 
nonzero temperature, even though in this paper we shall 
only use it at T = 0. At nonzero T, the warp factor / 
appearing in the metric (3.1 ) is given by (13. 21). When we 



later want to recover T = 0, we will simply set / = 1. 



We define H 



MN 



L 2 



IAIN 



and 



4 f Af If 



(2uq 2 -f) 



<r* ,J </V! H k 



Aq 2 f 



H, 



05 



l 05 ■ 



(3.38) 



where the zeroth and fifth coordinates are t and u respec- 
tively and where i and j run over the three spatial co- 
ordinates. The quantity Z clearly transforms as a scalar 
under rotations. With some effort, one can show that Z 
i s inva riant under the infinitesimal gauge transformations 
(3.37). Z is therefore a suitable gauge invariant quantity 



for our purposes. 



Five degrees of freedom in /ijvfJV can be eliminated by exploiting 
the gauge freedom to set h 5 M = for all M, reducing the num- 
ber of degrees of freedom in tiMN to ten. However, just as in 
electromagnetism, this does not completely fix the gauge. There 
exists a residual gauge freedom which allows one to eliminate five 
more components of h^jN on an y u = const, slice. 



There exist many different helicity zero, one and two gauge 
invariant combinations of /iMJVi but only five are independent. 
Different gauge invariants of the same helicity are related to each 
other by the linearized Einstein equations. 



12 



The dynamics of Z are governed by the linearized 



Einstein equations (3.31|. Using the linearized Einstein 



equations, it is straightforward but tedious to show that 
Z satisfies the equation of motion 



where 

A = 

B = 



AZ' + BZ = S . 



24 + V?i 2 + 6/ + q 2 u 2 f - 30f 2 

uf (u 2 q 2 + 6-6/) ' 
2 , q 2 u 2 (U-5f-q 2 u 2 ) + 18(4-/-3/ 2 ) 



P 



u 2 f (q 2 u 2 + 6-6/) 



S 8 , 4(q 2 u 2 +6-6/) 
— T^nn + 



/ 



00 



3uq 2 f 



(g 2 <J ij '-3gV>ij 



(3.39) 

(3.40) 
, (3.41) 
(3.42) 



+ 



, Su [g 2 (g 2 u 2 +6) - / (l2g 2 -9/'Q] 
/ 05 + " 3/ 2 (<? 2 U 2 - 6/ + 6) 00 



■^55 — &iq l ti5 . 



The connection between Z and the energy density may 
be found by considering the behavior of Z and Hmn 
near the boundary. Choosing for convenience the gauge 
H$m = 0, one can solve the linearized Einstein equations 



(3.31 ) with a power series expansion about u — in order 



to ascertain the asymptotic behavior of the metric per- 
turbation. Setting the boundary value of to vanish 
so the boundary geometry is flat and considering sources 
tjtfN corresponding to strings ending at u = 0, one finds 
an expansion of the form [25l [29] 



(3.43) 



H, v (u) = H$u 3 + H$u i + 



In the gauge H 5 m — the variation of the gravitational 
action (3.341 relates the asymptotic behavior of to 
the perturbation in the boundary energy density via [17] 



£ 



ji -"oo 

K 5 



(4) 



(3.44) 



The coefficient Hqo can m t urn be related to the asymp- 



totic behavior of Z by substituting the expansion (3.43) 



into Eq. ( |3.38[ ). In doing so one finds that Z has the 
2) u 2 + Z {3) u 3 + • ■ ■ , (3.45) 



asymptotic form 
Z(u) 
and that 



Z, 



H, 



(4) 
00 



16 



(3) 



(3.46) 



We therefore see that the energy density is given by 

; 5 



(3.47) 



The coefficient Z^) in the expansion (3.45), which has 
delta function support at the location ot the quark, gives 



the divergent stress of the infinitely massive test quark. 
It is therefore not of interest to us. We now see explicitly 
that, as we argued above, in order to obtain the energy 
density in the boundary quantum field theory the only 
aspect of the metric perturbation that we need to com- 
pute is Z , and furthermore that all we need to know are 
the coefficients in the expansion of Z about u = 0. 



D. The solution to the bulk to boundary problem 
and the boundary energy density 

Although we have defined Z at nonzero temperature, 
henceforth as we determine Z we return to T = 0, mean- 
ing / = 1. At zero temperature the coefficients A and B 
appearing in (3.39) are given by 



A = — 



B = U! 2 - q 2 



9 



(3.48) 
(3.49) 



and the general solution to (3.39) may be written 
u 3 I (uQ) 



Z(u) 



, K (u'Q) 

au — biu ) + a 

w 



u 6 K (uQ) lim 



du 



12 

U' 2 



S(u') + P 



(3.50) 



where Q = \J q 2 — uj 2 , Iq and Kq are modified Bessel 
functions, and a and (3 are constants of integration. 

The constants of integration are fixed by requiring that 
Z(u) satisfy appropriate boundary conditions. As Iq(uQ) 
diverges as u — > oo, regularity at u = oo requires a = 0. 
The constant j3 is fixed by the requirement that the e — > 
limit exists (so all points on the string contribute to 
the induced gravitational disturbance) and that no log- 
arithms appear in the expansion of Z(u) near u = 0. 
This last condition is equivalent to the boundary condi- 
tion that the metric perturbation Hmn vanish at u — 
so the boundary geometry is flat and unperturbed. For 
strings which end at u = 0, we have 



S(u) = s +0{u 2 ). 



where 



so 



lim u d, 

5 - >o \ u 



2 



a 2 v 00 
-OKi; lim — 

u-¥0 U 



(3.51) 



(3.52) 



where we have used the fact that too oc it at small u. 
Furthermore, the Bessel functions have the asymptotic 
expansions 



I a (uQ) = l + (D(u 2 ), 
K (uQ) = -7e - log(|uQ) 



0{u 2 ). 



(3.53) 
(3.54) 
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where 7e is Euler's constant. Substituting these expan- 
sions into the second integral in (3.50), we find 



,I {u'Q) , 
du — S(u ) 



Therefore, in order for the e 
have 



•S() 



Oil). 



(3.55) 



limit to exist we must 



so 
e 



(3.56) 



As one may readily verify, this value of /3 also eliminates 
all logarithms appearing in the expansion of Z near the 
boundary. 7 

With constants of integration determined and with the 
asymptotic forms (3.51 )— (3.54), it is easy to read off from 



find 



(3.50) the asymptotic behavior of Z. In particular, we 

(\ 

duQ(u)S(u) — sq I - + eQ(e) 



Zin\ = lim 
(i> e-s-0 



where 



(3.57) 
(3.58) 



is the bulk to boundary propagator. 

At this point it is convenient to Fourier transform back 
to real space. At zero temperature the source (3.42) 
Fourier transforms to 



S(t, x, u) 



s ,rX ( *™\ - ^V 2 (2/ ( 



3 



2^55 "f~ ta ) 

8VA5 , (3-59) 



where sums over repeated spatial indices i and j are 
implied. To take the Fourier transform of the bulk to 
boundary propagator one must give a prescription for in- 
tegrating around singularities at u) = q. The relevant 
prescription comes from causality and requires analyt- 
icity in the upper half frequency plane. This requires 
sending u> — > u> + ie. With this prescription we have 

G(t, r, u) = -^9(t)8'(-t 2 + r 2 + u 2 ) . (3.60) 

We therefore have 

//■oo 
dV J dug{t-t' \r-r' \u)S{t' >' » 

s (t,r) 



es (t',r')g(t',r-r',e) 



(3.61) 



7 A deformation in the boundary geometry implies that logarithms 
will occur in the expansion of Z at order n 3 and beyond. The lin- 
earized Einstein equations relate all higher order log coefficients 
to that of the u 3 coefficient in a linear manner. Therefore, if the 
order u 3 coefficient vanishes all other coefficients vanish. As is 
easily seen from the asymptotic forms ( | 3 ■ 5 1 [ ) - l |3.54[ l, the value 
of /3 given in l |3.56| l ensures the order u 6 logarithm coefficient 
vanishes. 



When we use (3.59) to express S(t',r',u) we find that 



when the du-integral of the first term in (3.59) is done 



by parts, the resulting boundary term cancels the e- 
dependent bulk-to-boundary propagator in (3.61). 



Upon assuming that the stress-energy tensor on the 
boundary does not have support at (t,r), i.e. that the 
observer is located away from the source, we can safely 
take the e — > limit and obtain for the energy density 



£(t,r) 



7T 



dV / du 0(t - t') 



(4ut 00 - iiwsVVW) 
+|r - rf (4*oo - 

- (*«v<wv;w) 



6"{W) 



At 55 + 2t H ) 
5"'(W) 



S"'(W) 



3 u 
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where 



(t-t') 2 + u 2 + \r-r'\ 



'|2 



(3.62) 



(3.63) 



-/Af 



and Vw is the partial derivative with r espect to X 
{£', r', u}. We see that the k| factors in (|3.47|) and (|3.59[) 



have cancelled, as expected. Because we have analytic ex- 
pressions for the profile of the rotating string, it turns out 
that we can evaluate £ explicitly, in two stages. First, us- 
ing (3.36) including its three-dimensional delta function, 



and using our expressions (3.17) and (3.19) for R(u) and 
4>(u), we obtain 



£(t,x) 



-n-2 



du 



with 



dt'[AS"(Ws)+BS'"(yv s )] , 

(3.64) 



A = 7 (3 + w 2 7 3 WqU (— t' + t + wy)) — u7Wo r sin cos ?/) 

+f 3 7 u! 2 ur sin 9 sin ip , 
B = « 4 7 V (1 - v 2 ^ 2 ^ 2 ) 

1 



+ I w 4 7 5 w 2 u 2 + 7 + 



37 



r s {t',u)Y 



-w 2 7(l + w 2 7 4 cj 2 u 2 )r 2 sin 2 -0 sin 2 9 
+2v 2 "/ 2 uj u r 2 cos ip sin ip sin 2 9 , 



2w 3 7 3 u;oM 2 (cos"0 — v 2 j 3 u)ousinip) rsir 

2„Vi 1 „.2.,4, ,2„.2\2 • 2 , 



where if) = ip + ^0(71* — t') and 

W s = -(t - t') 2 + u 2 + \r- r s (t', u)\ 2 



(3.65) 



(3.66) 



Next, the integral ( 3.64 ) can be carried out via the change 
of variables 



1 



1 



(7« -t'+t). (3.67) 



Since points with —(t'—t) < ju are causally disconnected 
from an observer on the boundary at time t, the argument 
of the (5-functions in Eq. (3.64) is non- vanishing. This 
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allows us to deform the domain of integration and the 
integral takes the generic form 



t — "fU 



dt'F n (t',u)S^(W s ) 
dv I d(F n ((,v)S( n \m), (3.68) 



In = du 

Jo 

2 r°° 

7 Jo 

where F n ((,i>) = F n (t'((, v), u((, v)). Furthermore, the 
argument of the ^-functions takes the form 

W s = r 2 + Rl + C (4f - 2vr sin 6» sin(2i/w + ip - oj Q t)) 
—2tRq sin0[^o;o sin(2z/u;o + <p — i^ot) 

+ cos(2t/o;o + <p — uj t)] , (3.69) 

which is then linear in £ and therefore allows us to eval- 
uate /„ formally, obtaining 



2 

2 

7 



riu 

(d c W) n+1 



^ 

Jo 



(-1)' 
i ( W 



d l c - l F n {C,v) 



x S^-V (W s (t' = t-2v,u = Q)). (3.70) 

Here Co{ v ) is the unique zero of W s at a given v and 
h> = (t — t rct )/2, but we will not describe in detail how 
these arise because the integrand in the first term in /„ 
vanishes: by direct inspection of F2 = A and F3 = B we 
find dlF 2 {C,v) — and d^,F 3 (C,iy) — 0. Conseq uently 
only the boundary terms in the second line of Eq. (3.70) 



contribute and we can perform the remaining integral 
using the (5-function. Finally, we obtain 



£ = 



2r 2 H 2 +4r7 2 H(t ret -t) 



247r 2 7 4 r 6 £ 6 
+ (2 7 2 -4rV 7 2 u; 2 sin 2 (9+3r 2 7 4 u; 2 S 2 )(i ret -t) 2 

+7n 2 u 2 E(t Ict -t) 3 + 4 7 2 w 2 (t rct -t) 4 
+8v"/ 2 tsj r(t Tet -t)(t rct -t+r S) 



x sin#cos((/9— uj t le t) 



(3.71) 



where S was given in ( 2.9 1. Recalling the discussion after 



(2.10), we note that £ is invariant under a shift of ip by 



(Hp that is compensated by shifts in both t and t rct by 
dip/ajQ. £ is a periodic function of <p with period 2tt and 
a periodic function of t with period 2tt/ujo, as it must be. 



Eq. (3.71) is the main result of our paper. It is an 
explicit analytic expression for the energy density of the 
radiation emitted by a test quark in circular motion, valid 
at all times and at all distances from the quark. The only 
aspect of the calculation that requires numerical evalua- 
tion is that one must solve the transcendental equation 



(2.10) for t Te t- With t YCt in hand, one uses (2.9) to com- 



pute c and then evaluates the energy density £ using the 
solution (13. 711). We shall illustrate our result 13.711) in 



several ways in Section IV. 



E. Far zone and angular distribution of power 

We close this section by evaluating the energy density 
( |3.71 ) in the r — > 00 limit and extracting the angular 
distribution of power a la Eq. (2.19). In the far zone 
limit we can replace S by its r — ¥ 00 limit, namely £ in 
(2.21 ), and we can safely replace t re t by (t—r) everywhere 



in (3.71 ) except within the (^-dependent argument of the 
cosine. In the far zone, the energy density (3.71) then 
reduces to that of Eq. (2.18) with 



j 2 V\ 4 - 4v 2 sin 2 9 - 7£ + 3£ 2 7 2 



24tt 2 



£ 6 7 2 



(3.72) 



distribution of power is given by 



Using Eq. (2.19), we find that the time-averaged angular 

(3.73) 



dP 

dn 



V 2 LU 2 ^f\ 



2 + v 2 sin 2 1 



16vr 2 



(1 



o • 2 /,\5/ 2 ' 

u 2 sm 2 6M 



which is the same as the weak-coupling result (2.24) up 
to a ^-independent overall factor! 

Upon integrating over all solid angles, we find the total 
power radiated 



P=^a 2 
2tt 



(3.74) 



where again a — u 7 2 ojo is the quark' s pr oper acceler- 
ation, reproducing Mikhailov's result (1.2) again. The 



fact that the power that we have obtained in this sec- 
tion by integrating over the angular distribution of the 
radiation (3.73) in the quantum field theory matches the 
energy flux (3.25) flowing down the classical string in the 



dual gravitational description is a nontrivial check of our 
calculations. 



IV. RESULTS AND DISCUSSION 
A. Radiation at strong coupling, illustrated 

Fig. [4] shows two cutaway plots of the energy density 
£ at strong coupling (multiplied by r 2 / P where P is the 
total power radiated) produced by a quark moving on a 
circle of radius Rq at velocities v = 1/2 and v = 3/4. 
The figure is obtained by evaluating (3.71). The motion 
of the quark is confined to the plane z = and at the 
time shown the quark is at x = Ro, y — 0, and is rotating 
counter clockwise. The cutaways in the plots show the 
energy density on the planes z = 0, <p = and ip = 7n/5 



where tp is the azimuthal angle. As is evident from the 
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FIG. 4: Left: a cutaway plot of r 2 £/P for v — 1/2. Right: a cutaway plot of r 2 £/P for v = 3/4. In both plots the quark is at 
x = i?o, y — at the time shown and its trajectory lies in the plane z — 0. The cutaways coincide with the planes z = 0, ip = 
and 93 = 77r/5. At both velocities the energy radiated by the quark is concentrated along a spiral structure which propagates 
radially outwards at the speed of light. The spiral is localized about — tt/2 with a characteristic width 80 ~ I/7. As i) — > 1 
the radial thickness A of the spirals rapidly decreases like A ~ 1 /j 3 . 
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FIG. 5: Plot of r 2 £/P at 9 = vr/2 and ip = 5tt/4 at t = as a function of r for u = 1/2. The plot illustrates the fact that 
the pulses of radiated energy do not broaden as they propagate outward. This implies that they do not broaden in azimuthal 
angle, either. Strongly coupled synchrotron radiation does not isotropize. 



figure, as the quark accelerates along its trajectory energy 
is radiated outwards in a spiral pattern. This radiation 
falls off like 1/r 2 and hence has a constant amplitude in 
the figure and propagates radially outwards at the speed 



of light. The figure shows the location of the energy 
density at one time; as a function of time, the entire 
pattern of energy density rotates with constant angular 
frequency Wq = RqV. This rotation of the pattern is 
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equivalent to propagation of the radiation outwards at 
the speed of light. As seen by an observer far away from 
the quark, the radiation appears as a short pulse just like 
a rotating lighthouse beam does to a ship at sea. 

What we see in Fig. [4] looks like an outward going 
pulse of radiation that does not broaden as it propagates. 
Analysis of (3.71 ) confirms this, as we have illustrated in 



Fig.[5]by extending the plot of the energy density outward 
to much larger radii at one value of the angular coordi- 
nates, at time t — 0. As time progresses, the pulses move 
outward at the speed of light, and an observer at large r 
sees repeating flashes of radiation. This figure provides a 
convincing answer to a central question that we set out to 
answer in Section I. We see that even though the gauge 
theory is nonabelian and strongly coupled, the narrow 
pulses of radiation propagate outward without any hint 
of broadening. We find the same result at much larger 
values of 7, where the pulses become even narrower - 
their widths are proportional to I/7 3 , as at weak cou- 
pling and as we shall discuss below. The behavior that 
we find is familiar from classical electrodynamics — in 
which the radiated energy is carried by noninteracting 
photons. Here, though, the energy is carried by fields 
that are strongly coupled to each other. And, yet, there 
is no sign of this strong coupling in the propagation of 
the radiation. No broadening. No isotropization. 

Fig. [5] is also of some interest from a holographic point 
of view. The qualitative idea behind gauge/gravity du- 
ality is that depth in the 5th dimension in the dual grav- 
itational description corresponds to length-scale in the 
quantum field theory. In many contexts, if one com- 
pares two classical strings in the gravitational descrip- 
tion which lie at different depths in the 5th dimension, 
the string which is closer to the boundary corresponds 
to a thinner tube of energy density in the quantum field 
theory while the string which is deeper, farther from the 
boundary, corresponds to a fatter tube of energy density. 
Our calculation shows that this intuitive way of thinking 
about gauge/gravity duality need not apply. The rotat- 
ing string falls deeper and deeper into the 5th dimension 
with each turn of its coils and yet the thickness of the 
spiral tube of energy density in the quantum field theory 
that this string describes changes not at all. 

The behavior of the outgoing pulse of radiation illus- 
trated in Fig. [5] is different than what one may have ex- 
pected for a nonabelian gauge theory given that solutions 
to the classical field equations in SU{2) and SU{2>) gauge 
theory are chaotic with positive Lyapunov exponents and 
are thought to be ergodic .'?() Of course, we have 

not solved classical field equations; we have done a fully 
quantum mechanical analysis of the radiation in a non- 
abelian gauge theory, in the limit of large N c and strong 
coupling. It is nevertheless surprising from the gauge 
theory perspective that the radiation we find turns out 
to behave (almost) like that in classical electrodynam- 
ics. From the gravitational perspective, we saw that the 
equations describing the five-dimensional metric pertur- 
bations linearize, meaning that there is no possibility of 
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FIG. 6: The energy density at 6 = tt/2 and ip — 5n/4 for 
quark velocity v — 1/2, as in Fig. [5] r ~ 33i?o corresponds to 
the location of a spiral in the energy density. Directly ahead 
of and directly behind the spiral, the energy density is slightly 
negative. To see how slightly, compare the vertical scale here 
with that in Fig. [5] 



chaotic or ergodic dynamics. 

The synchrotron radiation that we have found at 
strong coupling is, in at least one respect, qualitatively 
different from that in familiar classical electrodynamics. 
The difference is in principle visible in Fig. [5j but the ef- 
fect is small and hard to see without zooming in — which 
we do in Fig. [6] For large enough v we find that the en- 
ergy density is negative in some regions of space! Directly 
ahead of and behind the spiral, the energy density dips 
slightly below zero. This is impossible in classical electro- 
dynamics. The effect that we have found is numerically 
small, as a comparison of the vertical scales in Figs. [6] 
and [5] makes clear, but it is nevertheless in stark con- 
trast to the weak coupling results of Section [TTJ where 
the energy density is always positive, as in classical elec- 
trodynamics. This small effect serves as a reminder that 
the calculation that we have done is quantum mechani- 
cal. In a quantum field theory the energy density need 
not be positive everywhere — only its integral over all 
space is constrained to be positive. In a strongly coupled 
quantum field theory, quantum effects should be large. 
So, it is not surprising that we see a quantum mechani- 
cal effect. What is surprising is that it is so small, and 
that in other respects Figs. [4] and [5] look so similar to 
synchrotron radiation in classical electrodynamics. 



B. Synchrotron radiation at strong and weak 
coupling 

Given the qualitative similarities between Figs. [3] and [5] 
and the physics of synchrotron radiation in classical elec- 
trodynamics and weakly coupled Af = 4 SYM theory, 
we shall attempt several more quantitative comparisons. 
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First, as can be seen from Fig. [4j the radiation emitted 
by the quark is localized in polar angle about the equa- 
tor, 9 = tt/2. From the far-zone expression for the energy 
density, Eqs. (2.181 and (3.72), it is easy to see that the 



characteristic opening angle about 9 = tt/2 scales like 
89 r~j I/7 in the v — > 1 limit. Furthermore, from Fig. [4] it 
can be seen that the radial thickness A of each pulse of ra- 
diation decreases as v increases. Again, from Eqs. ( |2.18[ ) 
and (3.72) it is easy to see that A ~ I/7 3 in the v 1 
limit. 8 Intriguingly, the scaling of both 59 and A with 7 
are identical to those at weak coupling. 

Emboldened by the agreement between the 7-scaling 
of A and 89 at weak and strong coupling, we have com- 
pared the shape of the spiral of energy density in the 
z = plane, for example as depicted at two velocities in 
Fig. |4j with the shape of the classic synchrotron radiation 
spiral ( |2.12 ). We find precise agreement, as illustrated in 
Fig. [7] At strong coupling, the energy density £ of (3.71 ) 
is proportional to 1 /S 6 and so has maxima where 3 has 
minima. We have already seen that in the r — > 00 limit, 
E — > £ of (2.21) whose minima lie on the spiral (2.13), 
which is the large- r approximation to ( |2.12 ). It can also 
be shown that the minima of S lie on the spiral (2.12 ) at 
all r. 



Recall that our derivation of (2.12) in Section II was 



purely geometrical, relying only on the fact that the ra- 
diation is emitted tangentially, in the direction of the 
velocity vector of the quark and the fact that the pulse 
of radiated energy density propagates at the speed of 
light without spreading. We have seen in Fig. [5] that at 
strong coupling and in the limit of a large number of 
colors, the radiation emitted by a rotating quark does 
indeed propagate outwards at the speed of light, with- 
out spreading. This justifies the application of the geo- 
metrical arguments of Section II to the strongly coupled 
radiation. The agreement with (2.11) and (2.12) then 



implies that the strongly coupled radiation is also emit- 
ted in the direction of the velocity vector of the quark: 
if it were emitted in any other direction, A would have 
the same 7-scaling as a and 89; if the radiation were, 
for example, emitted perpendicular to the direction of 
motion of the quark, the spiral of energy density would 
have the shape ( |2.14 | instead of ( |2.12[ ) — see Fig. [7j We 
therefore reach the following conclusions: at both weak 
and strong coupling, the energy radiated by the rotating 
quark is beamed in a cone in the direction of the velocity 



8 There may be a holographic interpretation of A ~ I/7 3 . We 
saw in Section |III A| that there is one special point on the string, 
namely the worldsheet horizon at u = u c . Using the correspon- 
dence between u and length-scale in the boundary quantum field 
theory, we expect u c to translate into a length scale in the rest 
frame of the rotating quark, corresponding to a length scale 
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FIG. 7: r 2 £/P for v = 1/2 in the z = plane. The color code 
is the same as in Fig. [4| with zero energy density blue and 
maximal energy density red. The spiral curve marked with 
the black dots is ( |2.12[ ), namely the place where the spiral 
of synchrotron radiation would be in electrodynamics or in 
weakly coupled M = 4 SYM theory. We see that the spiral 
of radiation in the strongly coupled gauged theory with a 
gravitational dual is at the same location. This indicates that, 
as at weak coupling, strongly coupled synchrotron radiation 
is beamed in the direction of the motion of the quark. For 
reference, the solid black line is ( |2.14 1, namely the place where 
the synchrotron radiation would be if the quark were emitting 
a beam of radiation perpendicular to its direction of motion. 



of the quark with a characteristic opening angle a ~ I/7; 
at both weak and strong coupling, the synchrotron radia- 
tion propagates outward in a spiral with the shape ( |2.12[ ) , 
with pulses whose width A ~ 1 /j 3 does not broaden. 9 

We now turn to the time-averaged angular distribu- 
tion of power radiated through the sphere at infinity. 
By time-averaging, we eliminate all dependence on az- 
imuthal angle but nontrivial dependence on the polar an- 
gle 9 remains. Fig.[8]shows the normalized time-averaged 
angular distribution of power in classical electrodynam- 
ics (2.23a) and strongly coupled SYM (3.73) for velocity 



in the inertial frame in which the center of motion is at rest. It 
is tempting to identify <5 with A at v — y 1. 



v = 0.9. As is evident from the figure, in both classical 
electrodynamics and Af = 4 SYM the power is localized 



9 Note that if the radiation is isotropic in the instantaneous rest 
frame of the quark, then in the inertial "lab" frame it will be 
beamed in a cone with opening angle ~ I/7 pointed along the 
velocity of the quark. This, together with the result that the 
pulses of radiation pro pagat e without broadening, would yield a 
spiral with the shape 1 2.12 1. 
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Classical electrodynamics 
Strongly coupled SYM 




FIG. 8: The normalized time-averaged angular distribution 
of power radiated to large distances by a quark moving in 
a circle with velocity v = 0.9 in classical electrodynamics 
and in strongly coupled Af = 4 SYM theory. This normalized 
distribution is identical in weakly and strongly coupled Af = 4 
SYM theory. In classical electrodynamics and in both weakly 
and strongly coupled Af = 4 SYM theory, the time averaged 
angular distribution of power is localized about 8 = n/2 with 
a characteristic width 59 ~ I/7. 



about 6 = it/2 with a width 89 ~ 1/7. The angular 
distribution of power is slightly more peaked for strongly 
coupled Af = 4 SYM. 

Last, it is amazing to compare the time-averaged angu- 
lar dis tribution of power in strongly coupled Af = 4 SYM , 
(I3.73I), to that in weakly coupled TV = 4 SYM, 12.24b. 



The two results are identical, up to an overall normaliza- 
tion. Not only is there no broadening or isotropization, 
the pattern of radiation is identical at strong coupling to 
that at weak coupling! Just as for the total power radi- 



ated, the angular distribution at strong coupling (3.73) 
is related to that at weak coupling (2.24) by making the 



substitution (2.26) in their prefactor. We know that the 



small difference between the normalized time-averaged 
angular distribution of power in classical electrodynam- 
ics and in weakly coupled Af = 4 SYM theory seen in 
Fig. [8] comes from the fact that a quark in Af = 4 SYM 
radiates scalars as well as vectors. It seems evident that 
the same is true for radiation in strongly coupled Af = 4 
SYM theory as well. Furthermore, in any strongly cou- 
pled conformal quantum field theory with a classical dual 
gravity description, the time-averaged angular distribu- 
tion of power is giv en by our strongly coupled Af = 4 
SYM theory result (3.73). It is remarkable how quanti- 



tatively similar radiation is at weak and strong coupling. 



C. Relation with previous work 

Next, we return to the discussion of the relation be- 
tween our present work and previous works that we began 
in Section I. In Ref. [3], states produced by the decay of 



off-shell bosons in strongly coupled conformal field theo- 
ries were studied. It was shown in Ref. [9] that the event- 
averaged angular distribution of power radiated from the 
decay is isotropic (in the rest frame of the boson where 
the total momentum vanishes) and independent of the 
boson's spin. This should be contrasted with similar 
weak coupling calculations in QCD, where the boson's 
spin imprints an "antenna" pattern on the distribution 
of radiation at infinity [TUl412j . Our results are not in- 
consistent with those of Ref. [9] . First of all, we and the 
authors of Ref. [pj are considering very different initial 
states. Second, as discussed in Ref. [pj, if their initial off- 
shell boson has a non-trivial distribution of momentum, 
this distribution imprints itself on the angular distribu- 
tion of radiation at infinity and generically results in an 
anisotropic distribution of power. It is only the spin of 
the boson which doesn't produce any anisotropy. 

In Ref. [13] it was argued that isotropization must oc- 
cur in the strongly coupled limit due to parton branch- 
ing. Heuristically, the authors of Ref. [T3] argued that 
at strong coupling parton branching is not suppressed 
and that successive branchings can scramble any initially 
preferred direction in the radiation as it propagates out 
to spatial infinity. However, our results clearly demon- 
strate that this need not be the case. We see no evi- 
dence for isotropization of the radiation emitted by the 
rotating quark. In fact, the angular dependence of the 
time-averaged angular distribution of power is the same 
at weak and strong coupling. And, just as for weakly 
coupled radiation, the strongly coupled radiation that 
we have illustrated in Figs. [4] and [5] propagates outward 
in pulses that do not broaden: their extent in azimuthal 
angle at a given radius remains the same out to arbitrar- 
ily large radii and, correspondingly, their radial width is 
also unchanging as they propagate. It is not clear in any 
precise sense what parton branching would mean in our 
calculation if it occurred, since no partons are apparent. 
But, the essence of the effect is the transfer of power from 
shorter wavelength modes to longer wavelength modes. 
In our calculation, that would correspond to the pulses 
in Fig. [5] broadening as they propagate. This we do not 
see. 



D. A look ahead to nonzero temperature 

It is exciting to contemplate shining a tightly colli- 
mated beam of strongly coupled synchrotron radiation 
into a strongly coupled plasma. Although this will have 
its failings as a model of jet quenching in the strongly 
coupled quark-gluon plasma of QCD, so too do all extant 
strongly coupled approaches to this complex dynamical 
problem. The calculation that is required is simple to 
state: we must repeat the analysis of the radiation emit- 
ted by a rotating quark, but this time in the AdS black 
hole metric that describes the strongly coupled plasma 
that arises at any nonzero temperature in any confor- 
mal quantum field theory with a dual gravitational de- 
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scription. The calculation will be more involved, and 
likely more numerical, both because the temperature in- 
troduces a scale into the metric and because to date we 
have only a numerical description of the profile of the 
rotating string [19 . If we choose T -C 1/Ro and 7> 1 
then the criterion (1.7) will be satisfied and the total 



power radiated will be as if the quark were rotating in 
vacuum [TS] - Based upon our results, we expect to see 
several turns of the synchrotron radiation spiral before 
the radiation has spread to radii of order 1/T where it 
must begin to thcrmalize, converting into hydrodynamic 
excitations moving at the speed of sound, broaden, and 
dissipate due to the presence of the plasma. We will be 
able to investigate how the length scale or scales associ- 
ated with these processes vary with R and 7. Note that 
the typical transverse wavelengths of the quanta of radia- 
tion in the pulse will be of order Rq /j while their typical 
longitudinal wavelengths and inverse frequencies will be 
of order Rq/j 3 , meaning that we will have independent 
control of the transverse momentum and the frequency 
of the gluons in the beam of radiation whose quenching 
we will be observing. 

Since at nonzero temperature the coils of the rotat- 
ing string only extend outwards to some 7-dependent 
''max, it is natural to expect that this r max will cor- 
relate with at least one of the length-scales describing 
how the spiral of synchrotron radiation shines through 
the strongly coupled plasma. Our guess is that r max 
will prove to be related to the length scale at which 
the energy carried by the nonhydrodynamic spiral of 
synchrotron radiation as in vacuum converts into hy- 
drodynamic waves. If so, we can estimate the para- 
metric dependence of r max at large 7, as follows. For 
r 3> r max, the power P cx yf\a 2 cx ^/X^/Rq radiated 
by the rotating quark will be carried by long wavelength 
hydrodynamic modes whose energy density will fall off 
oc y/\^f 4 /r 2 . For r <C 1/T, we will have a spiral of syn- 
chrotron radiation whose energy density is e/r 2 with e 
given by (3.72). Recalling from Sect ion II that on the 
spiral £ = 1 — v cx I/7 2 , we see from (3.72) that, on the 



ated by a factor oc exp(— const rT). The parametric de- 
pendence of the r max at which the hydrodynamic modes 
take over from the nonhydrodynamic spiral can then be 
determined by comparing the parametric dependence of 
the two energy densities: 

~ \ JL exp (-const r max T) , (4.2) 



meaning that 

log 7 

^max OC ^ 

Interestingly, the numerical results of Ref . [TH] 
that r max diverges slowly as 7 — > 00. 
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spiral, the energy density is cx jr 1 . We expect that 
for 1/T -C r <C r max this energy density will be attenu- 
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